Equation of state and phase transitions in asymmetric nuclear matter 
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The structure of the 3-dimension pressure-temperature-asymmetry surface of equilibrium of the 
asymmetric nuclear matter is studied within the thermal Thomas-Fermi approximation. Special 
attention is paid to the difference of the asymmetry parameter between the boiling sheet and that 
of the condensation sheet of the surface of equilibrium. We derive the condition of existence of 
the regime of retrograde condensation at the boiling of the asymmetric nuclear matter. We have 
performed calculations of the caloric curves in the case of isobaric heating. We have shown the 
presence of the plateau region in caloric curves at the isobaric heating of the asymmetric nuclear 
matter. The shape of the caloric curve depends on the pressure and is sensitive to the value of 
the asymmetry parameter. We point out that the experimental value of the plateau temperature 
T ~ 7 MeV corresponds to the pressure P = 10 -2 MeV/fm 3 at the isobaric boiling. 



PACS numbers: 21.65.+f, 24. 10. Pa 



I. INTRODUCTION 



At zero temperature, T = 0, a saturated nuclear 
matter is stable at vanishing pressure P = 0. As the 
temperature increases, the equilibrium state is derived 
through the isotherm (equation of state) P(p,T), where 
p is the particle density. In general, there are two dif- 
ferent regimes of the equilibrium state, see Ref. @, Ch. 
9. If the evaporation time is large enough and the evap- 
orated particles are carried away from the heated liq- 
uid, the equilibrium condition is given by P(p,T) = 0, 
i.e., it has the same form as at T = 0. By increas- 
ing the temperature one can reach the so-called phase 
separation temperature T 8 (^-^] above which the equi- 
librium condition P(p, T) — is not satisfied anymore. 
This happens if the minimum of the isotherm P(p, T) is 
shifted on the p-axis and it is located at p — p mul ^ 0. 
The temperature T s is then the limiting one for the ex- 
istence of a stable state of the liquid phase under the 
equilibrium condition P(p, T) — 0. In another regime, 
the heated liquid is surrounded by the saturated vapor. 
The equilibrium condition requires then that the pres- 
sure, P hq (p,T), and the chemical potential, p hq (p, T), 
of the liquid phase should be equal to the corresponding 
ones, P va P(p, T) and p vap (p, T) , for the saturated vapor, 
see Ref. JjJ, Ch. 8. (Here and below indices "liq" and 
"vap" denote the liquid and vapor phases, respectively.) 
In this regime, as the temperature increases, the liquid 
becomes unstable with respect to the phase transition 
(boiling) and the corresponding boiling temperature is 
the limiting one for the heating of the nuclear liquid. 

The limiting temperatures and the critical behavior 
of the symmetric nuclear matter have been intensively 
investigated using an effective interaction and finite tem- 
perature selfconsistent mean- field theory The 
thermal properties of asymmetric nuclear matter have 
been considered in Refs. (| ||,|| |f|. In general, the 



two-component liquid and its saturated vapor co-exist 
with different asymmetry parameter X because of the 
p-dependence of the symmetry energy. The equilibrium 
state of the two phases of a two-component system is 
described by the 3-dimension surface of equilibrium in 
variables P,T and X, see Ref. @, Chs. 8 and 9. The 
evaluation of the (P, T, A)-surfaces of equilibrium and 
the analysis of the boiling regime for an asymmetric nu- 
clear matter is a main goal of the present paper. We 
also study the properties of the caloric curve in the two 
regimes described above. 



II. EQUATION OF STATE OF HOT NUCLEAR 
MATTER 

We will follow the temperature dependent Thomas- 
Fermi approximation using the Skyrme-type force as the 
effective nucleon-nucleon interaction. The energy den- 
sity, £ , and the entropy density, S, are given by H] 



+ 2*o[(l + V2V 2 - (so + l/2)(p£ + P 2 P )] 



+ 7^3P a [(l + x 3 /2)p 2 - (x 3 + 1/2) (pi + p% (1) 



S = (o A q J 3/2{Vq) ~ VqPq), 



(2) 



where x$, t% and a are the Skyrme force parameters, q is 
the isospin index (q = n for neutron and q = p for pro- 
ton), p q is the nucleoli density and p — p n +p p . The Fermi 
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integral Ju{Vg) = J dz z v /(l + exp(z — r) q )) in Eqs. (|l|) 
o 

and (Q) depends on the fugacity r\ q . The value of r\ q can 
be found from the condition 



Pq = A*Ji/ 2 (r} q ). 



(3) 



Here, A* = (1 / '2tt 2 ) (2m*T '/h 2 ) 3 / 2 and to* is the effective 
nucleon mass derived by 

h 2 H 2 1 
2 mi 2 rn 4 



+ 7 [t2 (a* + 1/2) - *i (xi + 1/2)] p q . 



(4) 



Using Eqs. (Q) and (g) one can derive the free en- 
ergy density J- = £—TS. Finally one obtains the pres- 
sure P (equation of state) and the chemical potentials 
fx q . Namely, 



P = P 2 



d_T_ 

dp p 



T,X 



dT 

dp n 



- ( d:F 
\dp P 



(5) 



where the asymmetry parameter X is defined as X — 
(p n — p p )/p. A numerical calculation of the pressure P 
with the energy density £ from Eq. ([I]) leads to the van 
der Waals-like isotherms P = P(T, X, pj, describing both 
the liquid and the vapor phases P,p|-|l(| . According to 
Gibbs's phase rule, the equilibrium states of the two- 
phases (liquid + vapor) and the two-components (neu- 
trons+protons) system are located on the 3-dimensional 
surface in the (P, T, X) space. To obtain the shape of 
the (P, T, A)-surface of equilibrium one should use the 
Gibbs equilibrium conditions jlj]: 



P li «(T,X,p)=P™?(T,X,p) , 



^{T,X,p) 



p v q ap (T,X,p) 



(6) 



We will imply that the nuclear matter is a bound sys- 
tem of nucleons assuming that the following additional 
condition p q lq < is also satisfied. 

III. SURFACES OF EQUILIBRIUM 

We have performed calculations of the surfaces of equi- 
librium in (P, T, X) space using the SkM force Jl4| with 
the parameters t = -2645 MeV-fm 3 , t x = 385 MeV-fm 5 , 
t 2 = -120 MeV-fm 5 , t 3 = 15595 MeV-fm 3+3<T , x = 0.09, 
x\ = X2 = X3 = and a = 1/6, adopted for the ground 
state of nuclei. We point out that the use of the Skyrme 



force from the ground-state calculations is a quite rea- 
sonable approximation for our purposes. The effective 
interaction is modified only by a few per cent in a wide 
temperature region from to 20 MeV (l^Jl^]. The nu- 
merical procedure we used starts from certain values of 
A lic i, T and looks for p U( i, p va ? and A va ? which satisfy 
the equilibrium conditions (^|). 
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FIG. 1. Surface of equilibrium in P, T, X-space. The up- 
per sheet is the surface of boiling and the lower one is the 
surface of condensation. The azeotropic line and the critical 
line are marked by the letters A and C, respectively. The 
drip line fi„ = is marked by the letter B. The point D indi- 
cates the maximum possible asymmetry for the bound liquid 
phase at \x n < and the point E is the critical point for the 
symmetric system X = 0. 

The section of the equilibrium surface calculated by 
our code is shown in Fig. [l]. This figure represents the 
equilibrium states having X > 0, i.e. states of the 
neutron-rich nuclear matter. We point out that there are 
two sheets of the surface of equilibrium. The upper sheet 
is the surface of boiling and the lower one is the surface 
of condensation. The interior space between sheets is the 
phase separation region where the two phases co-exist. 
The crossing points of both sheets with a straight line 
P = const, T = const, give the equilibrium asymmetry 
parameters A hq and A vap for liquid and vapor phases, 
respectively. In general the asymmetry parameters A hq 
and A vap are different. The vapor asymmetry A vap ex- 
ceeds the corresponding liquid one A llq (at A llq > 0). 
This is a feature of the nuclear matter. The density de- 
pendence of the isospin symmetry energy (see Eq. (0)) 
provides the condition \p n \ < \p p \ and the preferable 
emission of neutrons. 

The above mentioned sheets of the surface of equilib- 
rium coincide along the azeotropic (P, T)-line in X = 
plane denoted by letter A in Fig. |l[ The azeotropic line 



is cut off at the critical temperature T ( 



(sym) 



14.61 MeV 



(point E) of the symmetric nuclear matter derived by the 
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condition M 



dp / 3 A n 



d 2 P 
dp 2 



= 



(7) 



T,X=0 



Both sheets of the surface of equilibrium also coincide 
along the critical line marked by letter C in Fig. [I]. The 
critical line of an asymmetric nuclear matter is derived 
by the condition, see Ref. @, Ch. 9, 



dX 



T,P 



\dX 2 



(8) 



We point out that the critical temperature T cr derived 
by Eq. (§) at X ^ is different from the one, T c ( r sym) , 
obtained from Eq. (Q) . If one goes along the critical line 
C, the critical temperature T cr decreases and the corre- 
sponding pressure increases as the value of X increases. 

The line B in Fig. |l| is the line with p, n = 0. The cross- 
ing point of the critical line C with the line B (point D 
in Fig. [l]) provides the maximum possible asymmetry for 
the bound (p}* q < 0) liquid phase. This point is located 
at X = 0.68, T = 10.4 MeV and P = 0.26 MeV/fm 3 for 
the SkM interaction. Note that the cold nuclear matter 
at T = is bound at X < 0.31. Thus, the hot nuclear 
matter can exist (in the bound state) at higher asymme- 
try than the cold one. This feature of hot nuclear matter 
appears because of the increasing of the symmetry energy 
with temperature. 



IV. BOILING AND CALORIC CURVES 



Let us consider the isobaric (T, X) phase diagrams ob- 
tained as the cut of the equilibrium surface of Fig. |l| by 
plane P — P cx t = const. These diagrams are very use- 
ful for obtaining an insight into the liquid-vapor phase 
transition in hot asymmetric nuclear matter. The shape 
of (T, A")-diagram depends on the value of the pressure 
P cx t- The (T, V)-diagram contains the critical point T cr 
if the value of pressure P ex t exceeds the critical (maximal 

allowed) pressure pjr ym ' ) on the azeotropic line in Fig. [j]. 
For the SkM force used here, the value of the critical pres- 
sure obtained from (0) is P c ( r ym) =0.2109 MeV/fm 3 . In 

the opposite case P ex t < Pcr ym \ the (T, V)-diagram con- 
tains the point of equal concentration X = 0. In Figs. || 
and H we have plotted the (T, V)-diagrams at the pres- 
sures P cxt = 0.15 MeV/fm 3 and P oxt = 0.25 MeV/fm 3 , 
respectively. In each diagram we have two lines. The 
lower line (on the left) is the line of boiling and the up- 
per line (on the right) is the line of condensation. In 
Fig. 0, these lines meet at the point of equal concentra- 
tion (X = 0) and in Fig. ||, they meet at the critical 
point T cr . The space between the boiling and condensa- 
tion lines corresponds to the states of co-existing phases. 




FIG. 2. A cut of the (P, T, X)-surface of equilibrium by 



(sym) 



the plane P = 0.15 MeV/fm 13 < P< 



Let us consider the boiling of the asymmetric nuclear 
matter at the fixed external pressure P ex t see 
Fig. ||. With the increase of temperature T at a cer- 
tain asymmetry parameter X = A"^ q , the asymmetric 
nuclear matter begins to boil at temperature Ta on the 
line of boiling (point A in Fig. ||) . The boiling process is 
accompanied by the preferable evaporation of neutrons 
because of the condition < \p p \. The liquid phase 
is then shifted to a more symmetric state and the boil- 
ing temperature increases due to the evaporation of the 
less-bound particles. The final state of system depends 
on the regime of heating. In regime I we assume that 
the liquid is evaporated in a free space and the vapor is 
taken away, the system moves along the surface of equi- 
librium to the point F, i.e., reaches the symmetric state 



X = X 



liq _ 



at T = T n - 



< Tcr yn ^ on the azeotropic 



line. The boiling is finished then at fixed T 
X = 0. The asymmetry of the vapor, V vap , decreases 
with increasing T and disappears at T — T max - This 
means, for instance, that the yield of symmetric clus- 
ters at the condensation of the vapor phase has to grow 
with T. In regime II we assume that the vapor remains 
near the liquid (at fixed pressure P = P ox t, this means 
an evaporation into a closed but not a fixed volume) the 
system moves along the boiling line and reaches the point 
D in Fig. ||. The point D corresponds to a fully evapo- 
rated liquid at asymmetry X™ p of the vapor (point C on 
the surface of condensation) that equals to the starting 
asymmetry X^ of the liquid phase. We point out that, 
in general, the trajectory of motion of the liquid phase 
along the surface of equilibrium can be located beyond 
the plane P = const and Fig. || should be considered as a 
formal illustration of the behavior of the heating system 
for the regime II where both points A and D are taken at 
different pressures. An analogous analysis of the boiling 
of the asymmetric nuclear matter can also be done in the 
case of P cxt > Pir ym ^ , see Fig. 0. However the presence of 
the critical point on the (T, X) phase diagrams in Fig. || 
leads to a very specific effect of the retrograde condensa- 
tion, see also Refs. PlJl^]- If one goes along straight line 
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AB in Fig. || (at closed volume), the liquid starts to boil 
at temperature An increase of the temperature leads 
to an increase of the evaporation. However the evapora- 
tion begins to decrease at certain temperature T < Tb 
and the vapor disappears at temperature Tb < To 
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FIG. 3. The same as in Fig. g for P = 0.25 MeV/fm 3 
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The (T, X) phase diagrams allow us to study the shape 
of the caloric curve for the case of isobaric heating. We 
will consider the case of the evaporation in a closed (but 

not fixed) volume at the fixed pressure P < Pcr ym \ i.e., 
regime II in Fig. g. Let us introduce the volume fractions 
A hq and A vap of the liquid and vapor phases defined by 
A iiq = y-iiq/y an d A vap = V vap /V, where V liq and V vap 
are the volumes of the liquid and vapor phases, respec- 
tively, and V = V hq + V vap . The excitation energy per 
particle, E*/A, is given by 



E* 
~A 



A liq gHq(pHq ! X liq , T) + A vap £ vap (/9 vap , X vap , T) 



A liq pliq _|_ A vap ^vap 



£ liq (p hq ,X liq ,T) 



oliq 



(9) 



T=0 



where £ liq (p liq , X liq , T) and £ vap (p vap , V vap , T) are, re- 
spectively, the energy densities of the liquid and vapor 
phases, derived by Eq. (Q) and taken at the correspond- 
ing values of the particle density p and the asymmetry 
parameter X. Let us consider a certain asymmetry pa- 
rameter X — Xq. To calculate the excitation energy 
E*/A we note that A liq = 1, A vap = 0, and X Uq = A 
if T is below the boiling line and A liq = 0, A vap = 1 and 
X vap = Xo if T is above the condensation line in the 
(T, X) phase diagram of Fig. ||. In the case of co-existing 
phases (space between lines of the boiling and the con- 
densation in Fig. ||) the required values of A llq , A vap , p llq , 
p vap , X llq and A vap can be found using the following 
procedure. At given Xq 7 we go along the straight line 
T = const in both directions and get the cross points 
with the boiling and condensation lines which provide the 
asymmetry parameters X hq and A vap and the densities 



p hq and p vap for both phases. The volume fractions A llq 
and A vap are then determined by the following equations 



X n = 



^liq^liqj^liq _|_ ^vap^vap j^vap 
A liq pliq _|_ ^vap ^vap 



A liq + A vap = 1 



(10) 

Caloric curves determined by Eq. (||) for X = 0.2 
and pressures P = 10~ 3 MeV/fm 3 , 10~ 2 MeV/fm 3 
and 10 _1 MeV/fm 3 are presented in Fig. |[ The solid 
line at low values of the excitation energy E*/A cor- 
responds to the heating of the degenerate Fermi liquid 
with E*/A ~ T 2 . The solid line at high excitation en- 
ergy E* I A describes the classical Boltzmann's gas with 
E* I A = (3/2)T. The region of two phase co-existence is 
displayed by the dashed line. The caloric curve is a con- 
tinuous function of E* /A and has a break in its derivative 
at two points connecting the two-phase region with the 
corresponding single-phase regions. As seen from Fig. 
the plateau region corresponds to the two-phase region. 
A small increase of the dashed line with T is due to the 
motion along the boiling path AD in Fig. g. The value 
of the plateau temperature increases with the increase of 
pressure. Experimental observations show a nearly flat 
caloric curve with a temperature of about 7 MeV [[17| ■ If 
one could assume the process of isobaric heating for the 
description of the experimental data, the order of mag- 
nitude of the pressure should be 10 -2 MeV/fm 3 for this 
process. To obtain more accurate estimation an analysis 
of the contribution of both the surface and the Coulomb 
forces into the equilibrium condition (p) is needed P^,[18| . 
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FIG. 4. Caloric curves for the isobaric heating of the 
asymmetric nuclear matter in regime II for different pres- 
sure P shown near the curves at fixed asymmetry parameter 
X = 0.2. 

The asymmetry dependence of the shape of the caloric 
curve is displayed in Fig. ||by plotting two curves at X = 
0.1 and 0.3. The figure shows that the plateau tempera- 
ture is slightly sensitive to the asymmetry parameter. At 
low asymmetry the two-phase region of the caloric curve 
is flatter and it is shifted to the lower values of E* /A as 
compared to the case of high asymmetry 
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FIG. 5. The same as in Fig. ^| for the fixed pressure 
P = 0.01 MeV/fm 3 and different values of the asymmetry 
parameter X shown near the curves. 

Figure || shows results of the calculation of the caloric 
curve assuming the regime I, i.e. the case when vapor is 
taken away from the liquid during the boiling. The cal- 
culation was carried out for the values of pressure 10~ 3 , 
10~ 2 and 1CP 1 MeV/fm 3 . The solid lines correspond to 
the heating of the liquid from T = to T = Tx with 
asymmetry parameter X = 0.2. The regions where sat- 
urated vapor is present (but taken away) are displayed 
by dashed lines. It is seen from Fig. ^ that the value of 
the external pressure does not change the shape of the 
caloric curve calculated for regime I. Only the interval of 
E* I A is sensitive to the value of the pressure. 
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FIG. 6. Caloric curves for the isobaric heating in the case 
when vapor is being taken away from liquid (regime I). Cal- 
culations presented for values of pressure 10 -3 MeV/fm 3 (a), 
10" 2 MeV/fm 3 (b) and 10" 1 MeV/fm 3 (c). 



V. SUMMARY 

Starting from the temperature dependent Thomas- 
Fermi approximation with the effective SkM force, we 



have evaluated the 3-dimensional (P,T,X) surface of 
equilibrium of the asymmetric nuclear matter. The 
(P, T, X) surface of equilibrium contains two sheets which 
correspond to the boiling and the condensation. Both 
sheets coincide along the azeotropic line and along the 
line of the critical pressure marked by C in Fig. [y. In 
general the vapor asymmetry A vap on the surface of con- 
densation exceeds significantly the corresponding liquid 
asymmetry X hq (at X hq > 0) lying on the surface of 
boiling. This is due to the fact that the isospin symme- 
try energy causes a preferable emission of neutrons. The 
surface of equilibrium is restricted by the drip line (line 
B in Fig. |l|) where one has fj, n = 0. The (P, T, X) surface 
of equilibrium in Fig. |l| shows also that the hot nuclear 
matter exists (in a bound state) at higher asymmetry 
than the cold one due to the increase of the symmetry 
energy with temperature. 

Our analysis of the isobaric (T, X) phase diagrams (see 
Figs. and ||) shows that the process of isobaric boiling 
of the asymmetric nuclear matter depends on the value 
of the fixed pressure P ox t. In the case of low enough 
pressure P oxt < P c ( r sym) =0.2109 MeV/fm 3 the (T,X)- 
diagram contains the point of equal concentration X = 
and the process of boiling is accompanied by a decrease 
of the vapor asymmetry with increasing temperature T. 
In the case of high pressure P ox t > Per™ j the process 
of boiling leads to the very specific effect of a retrograde 
condensation with growingtemperature T, see Fig. |3| 

We have shown in Figs. Q and || that the caloric curve 
for the process of isobaric heating of the asymmetric nu- 
clear matter contains a plateau region where both the liq- 
uid and the saturated vapor phases co-exist. The position 
of the plateau region on the caloric curve depends on the 
pressure P cxt and is almost insensitive to a change in the 
asymmetry parameter X . As can be seen from Fig. ||, the 
shape of the caloric curve is significantly changed with 
asymmetry. At low asymmetry the two-phase region of 
the caloric curve is flatter. We pointed out that the ex- 
perimental observation of the saturation of the caloric 
curve at temperature of about 7 MeV is obtained in our 
approach as an isobaric boiling of the nuclear liquid at 
the pressure P ox t ss 10 -2 MeV/fm 3 . Note, however, that 
the effects of nuclear surface and Coulomb interaction 
are not included in our study. The asymmetry parame- 
ter X vap of the vapor phase decreases along the boiling 
path of the caloric curve (dashed line in Figs. || and ||). 
Thus, the present model predicts an increase of the yield 
of the symmetric clusters at the condensation of the va- 
por phase with an increase of the excitation energy E* / A 
in the plateau region of the caloric curve. 
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